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Abstract 

The tunneling approach to the wave function of the universe has been recently criticized 
by Bousso and Hawking who claim that it predicts a catastrophic instability of de Sitter 
space with respect to pair production of black holes. We show that this claim is unfounded. 
First, we argue that different horizon size regions in de Sitter space cannot be treated as 
independently created, as they contend. And second, the WKB tunneling wave function is 
not simply the 'inverse' of the Hartle-Hawking one, except in very special cases. Applied 
to the related problem of pair production of massive particles, we argue that the tunneling 
wave function leads to a small constant production rate, and not to a catastrophy as Bousso 
and Hawking's argument would suggest. 



1. Introduction 

The tunneling proposal for the wave function of the universe has been recently criticized 
by Bousso and Hawking |l|] who argue that it predicts a catastrophic instability of de Sitter 
space with respect to production of black hole pairs. Here it will be shown that their 
argument is incorrect. 

Let us first review the Bousso-Hawking argument. First, they estimate the pair produc- 
tion rate for black holes in an inflating universe as 

r oc expf-^B/^ + 5*0], (1) 

where 5*0 and Sbh are Euclidean actions for de Sitter and (degenerate) Sczwarzschild-de Sit- 
ter instantons, respectively. The Schwarzschild-de Sitter instanton, which has the topology 
of S2 X S2, is interpreted as describing nucleation of two black holes whose horizon radii 
are equal to the radius of the de Sitter horizon, H~^. A simple calculation gives (in Planck 
units) 

r ocexp(-7r/3i7^). (2) 

For H 1, black hole nucleation is exponentially suppressed. Evaluation of the tunneling 
rates using instanton techniques is a standard practice in the literature, and results similar 
to (|1]), (0) for black hole nucleation have been previously derived by Mellor and Moss 0] and 
others [Q. We have no objections against these results 

Now, Bousso and Hawking suggest that instead of using the standard instanton approach, 
the nucleation rate can be calculated directly from the wave function of the universe. They 
argue that each horizon volume in an inflating universe can be regarded as having been 
nucleated independently of other horizon regions. The rate of black hole production can 
then be found as 

F (X Pbh/Po, (3) 

where Pbh is the probability of quantum nucleation for a horizon-size universe containing a 
pair of black holes, and Pq is the same without black holes. 

The nucleation probability in quantum cosmology is sensitive to the choice of boundary 
conditions for the wave function of the universe. With Hartle-Hawking boundary conditions. 
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the wave function is given by the integral 

i^HH = J e-^ (4) 

where 5* is the Euchdean action and the integration is taken over compact Euchdean geome- 
tries and matter fields with a specified field configuration at the boundary. The boundary 
configuration is chosen to be the 3-geometry at the moment of nucleation. For a de Sitter 
universe this is the maximal 5*3, while for black holes in de Sitter it is 5*1 x 5*2. The proba- 
bility is found from P oc Assuming that the severe divergence problems of the integral 
(^) are somehow resolved, one can expect that the dominant contribution to is given by 
the instantons, 

ipHH OC e"^, (5) 

so that P oc e^^'^, where 25 is the corresponding instanton action. Then Eq. @ reduces 
to (|l]) and (@). This completes the quantum-cosmological derivation of the black hole pair 
production rate. 

Turning to the tunneling wave function, Bousso and Hawking argue that, since it is 
suppressed rather than amplified under the barrier, this wave function should be given by 

"Vt" oc e+^. (6) 

Then the sign of the exponent in (|l]) is changed to the opposite, and Eq. (§) is replaced 
by r oc exp(+7r/3if^). This shows no suppression of black hole production, indicating an 
instability of de Sitter space. 

We have two objections to this argument. First, we disagree that horizon volumes in an 
infiating universe can be treated as nucleating independently. Second, the tunneling wave 
function is not given by Eq. @, except in some very special cases. 

In the next Section we explain why we believe that different horizon volumes cannot 
be regarded as nucleating independently. In Section 3 we discuss the tunneling boundary 
condition for ip and explain why it is very difficult to implement in the problem of black 
hole nucleation. Fortunately, the tunneling wave function can be found in a very similar 
problem of massive particle production during inflation. The Bousso-Hawking argument can 
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be applied to this problem without change. If the argument were true, the tunneling wave 
function would predict a catastrophic instability of de Sitter space with respect to particle 
production. However, we show in Section 4 that it actually predicts a small, constant pair 
production rate, in agreement with instanton calculations. 

2. Black hole nucleation vs. nucleation of universes with black holes 

To see why different horizon volumes cannot be regarded as nucleated independently, 
consider a generic case of inflation driven by a scalar field with a slowly-varying potential 
V{(f)), disregarding for the time being the process of black hole production. The field 
is approximately constant on the horizon scale, and nucleation of horizon-size universes 
with different values of can be approximately described by de Sitter instantons. The 
corresponding probability is 



where H'^{4>) = 87rV(0)/3 and the upper and lower signs are for the Hartle-Hawking and 



functions do differ mainly by the sign in the exponent). Suppose an observer occupies a 
horizon-size region with = 0o and wants to know the probability distribution for in 
the adjacent regions. If different horizon volumes were nucleated independently, then this 
distribution would be given by (|^) and would be independent of 0o. This would imply rapid 
variation of the field from one horizon region to another. Of course, the observer will 
not be able to see this variation as long as infiation continues. But after infiation, all these 
horizon volumes will eventually come within the observer's horizon, and rapid variation of 
will manifest itself in large density fiuctuations on all scales. This conclusion is in a sharp 
disagreement with the standard analysis of density fiuctuations in the infiationary universe 



The evolution of the field during infiation is determined by the quantum state of the 
field. It is usually assumed that this state is close to the de Sitter-invariant (Bunch-Davies) 
vacuum. Then, it has been shown that the evolution of can be pictured as a diffusion 
process, resulting in strong correlations between the values of in adjacent horizon regions. 
As a consequence, density fluctuations are small on sufficiently small scales. It is possible 




(7) 



tunneling wave functions, respectively 0; 3- (This is one example when the two wave 
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that some highly excited state of the field would result in evolution suggested by the random 
nucleation picture. Hence, the quantum-cosmological problem of finding the probability for 
the universe to nucleate with a given value of is not the same as the problem of finding 
the probability distribution for within the universe. Quantum cosmology can be said to 
determine the latter distribution only in the sense that it determines the initial quantum 
state of at the nucleation of the universe. 

Quite similarly, the probability of black hole pair production in an infiating universe is 
determined by the quantum state of the gravitational field. It is not related in any simple 
way to the probability for the universe to nucleate with a pair of black holes. Eq. (^) 
with the Hartle-Hawking wave function ipHH correctly gives the leading exponent of the pair 
production rate F only because ipHH is defined in terms of the same Euclidean path integrals 
used in the instanton evaluation of F. 

3. The tunneling wave function 

The wave function of the universe ip is defined on superspace, which is the space of all 
3-geometries (and matter field configurations). It satisfies the Wheeler-DeWitt equation 



where 7i is a differential operator on superspace. The tunneling boundary condition, as 
defined in Refs. [||, |ll|], requires that ip includes only outgoing waves at singular boundaries 
of superspace. The physical meaning of this condition is that the universe originates in a non- 
singular way (at the regular boundary), but may end at a singularity. The regular boundary 
of superspace includes all singular 3-geometries which can be obtained as non-degenerate 
slices of smooth Euclidean 4-geometries. 

The word "non-degenerate" means that the presence and the nature of the singularity are 
stable with respect to small perturbations of the slicing. To give an example of a degenerate 
slicing, consider a torus lying on a horizontal surface and imagine slicing it with horizontal 
planes. The slice at the bottom is a circle. But if the torus is slightly tilted, the circular slice 
disappears, and the bottom slice is an isolated point. The circular slice is degenerate in the 
sense that it is present only for a very special slicing. A rigorous definition of non-degenerate 
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4 



slices is given in Morse theory |T^, |13| . It can be shown that singularities on non-degenerate 
slices always occur at isolated points. 

In addition to the tunneling boundary condition, one has to impose some kind of regu- 
larity condition, e.g., \ip\ < oo. 

We do not know how to solve the Wheeler-DeWitt equation (|^), except in simple models of 
high symmetry, when the infinite number of degrees of freedom of the universe can be reduced 
to one or two, with the remaining degrees of freedom being treated as small perturbations. 
For example, the nucleation of an inflating universe (without black holes) can be described 
by treating the radius of the universe, a, as a non-perturbative variable and deviations from 
spherical symmetry as perturbations. The regular boundary of superspace is then at a = 0, 
with finite amplitudes of the perturbations. 

To find the tunneling wave function for a universe with a pair of black holes, one needs to 
construct an appropriate minisuperspace model. The 5*2 x 5*2 instanton and the symmetry 
of the problem suggest a minisuperspace consisting of 3-geometries having topology Si x 5*2 
with the radius of the sphere 5*2 fixed at H~^. This model would have the radius r of the 
circle 5*1 as its only variable and would certainly be simple enough to solve. However, the 
minisuperspace boundary at r = cannot be regarded as regular boundary. This is clear 
from the fact that this singular geometry has singularities on a 2-dimensional surface (S'2), 
rather than at an isolated point . A more complicated model, where the radii of Si and 5*2 



are both allowed to vary, has a similar problem. The tunneling boundary condition cannot 
be implemented in such models. Moreover, it appears that such models are not suitable as 
minisuperspaces, since they include configurations which are unstable with respect to small 
perturbations. Extensions of Si x S'2 models that would not have this problem are bound 
to be less symmetric and far more difficult to solve. We note that even if we managed to 
construct a suitable minisuperspace model and solve for the tunneling wave function ipT, 
this result would be relevant only for the problem of the nucleation of a universe with a pair 
of black holes, and not for the rate of black hole production during inflation. 

4. Nucleation of massive particles 

Instead of trying to fix the Si x S'2 minisuperspace for black hole nucleation, we shall 
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consider a related problem of nucleation of massive particles. The corresponding instanton is 
the de Sitter S4 with a particle world line in the form of a big circle |T^. A similar instanton 
for nucleation of small extremal black holes was discussed by Mellor and Moss 0. The 
instanton action is Sp = Sq + 27rH~^m, and the nucleation rate is 

r oc exp[— + So] = exp(— 27rm/i7). (9) 

Particle nucleation in deSitter space has all relevant features of black hole nucleation, and 
the Bousso-Hawking argument can be applied to it without change. If the argument were 
true, then the tunneling wave function would predict a catastrophic instability of de Sitter 
space with respect to pair production. However, we will show that in fact it gives the correct 
particle production rate (H). 

We shall assume that the mass of the created particles is sufficiently small that their 
gravitational backreaction can be ignored, m << 1, yet sufficiently large that the concept 
of particle can be unambiguosly defined in our expanding background. As we shall see, this 
requires m » H. This is actually the same mass range in which the world line instanton 
approximation is valid. 

Within this regime, we can represent the particles as excitations of a massive scalar field 
(p. Let us expand the field configuration as (p = J2 fnQni where Qn are the harmonics on the 
3-sphere. Then the Wheeler-De Witt equation (H) can be solved perturbatively in using 
the WKB ansatz tfj ~ exp iS, where 

S = So{a) + \Y.Sn{a)fl (10) 

and 

^ = ±a{a'H^ - (11) 
da 

The solution of the Hamilton- Jacobi equation (|ll]) is the action along the classical de Sitter 
solution, which represents an infiationary universe. For a < classical motion is forbid- 
den, Sq is imaginary, and ipisa. linear combination of growing and decaying exponentials. For 
a > classical motion is allowed and is a linear combination of incoming and outgoing 
waves, whose fiux points towards and away from the forbidden region respectively. The tun- 
neling boundary condition at a — > 00 requires that only the outgoing wave should be present. 
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The matching conditions for ip then imply that the growing and decaying components have 
comparable magnitude near the classical turning point a = H~^. 

With our ansatz, the Wheeler-De Witt equation reduces to the functional Schrodinger 
equation for the quantum field in a fixed de Sitter background, with the scale factor 
playing the role of time. Its solution is given by Sn = H^a^iyn/^n, where a dot indicates 
derivative with respect to the conformal time r, defined by a = {H cost)~^ , and z/„(r) are 
normal modes of the classical scalar field equation. For each n there are two independent 
solutions, and the quantum state is specified once we choose a particular linear combination 
of these as our 'negative frequency' mode z/„. The mode functions Un should in principle be 
determined by the outgoing-wave and regularity conditions at the boundaries of superspace: 
a — > 0, fn—>- ±C)0. This, however, cannot be implemented within the perturbative approach, 
since the expansion ([T0|) breaks down at large values of /„. Instead, we shall follow Ref. [|1^ 



and require that the wave function does not increase towards large /„. It appears that this is 
the best one can do to represent the boundary condition i'dfnl — ^ oo) < oo. Mathematically, 
our condition requires that Im{Sn) > along the three branches of the semiclassical wave 
function (growing, decaying and outgoing). It can be shown that this requirement uniquely 
determines z/„, and hence the quantum state of the scalar field [|T^. The corresponding 
eigenmodes are given by 

Un oc cosre*"^F(/i, 1 -/i,n + 1; (1 + 2tanr)/2), (12) 
_ 1 _ /9 _m2y/' 

^ ~ 2 1,4 ' 

where r ranges from to 7r/2 for the outgoing branch and from to ±ioo for the growing and 
decaying exponentials. These modes correspond to the Bunch-Davies or de Sitter invariant 
vacuum. 

The Hartle-Hawking wave function for this model includes only the growing exponential 
for a < and a linear combination of incoming and outgoing waves with equal amplitudes 
at a > H~^. The mode functions z/„ are given by Eq. (0) and are the same as for the 
tunneling wave function. It is easily seen that the two wave functions are not related by the 
transformation S — > —S, as suggested by Bousso and Hawking [see Eqs. (|^), (|)]. 
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It is well known that a particle detector responds in the Bunch-Davies vacuum as if 
there was a thermal distribution of particles with temperature H/2tt [0. This is already in 
qualitative agreement with the instanton result (^. The correspondence can be made even 
more precise by using the method of Bogoliubov coefficients pO| . 

For m » H and for each n an unambiguous definition of particles can be given at late 
times. The scalar wave equation is given by 

i>n + 3if tanh(m) z>„ + rn?Vn + ^z^n = 0, (13) 

where now the dots indicate derivative with respect to cosmological time defined by a = 
cosh(i?t). For t >> t„ where t„ is the time at which the particles become non-relativistic, 
a(^n) ~ nH/m, the last term can be ignored and we have approximate solutions of the form 

z/„oca-3/2(a„e*-* + /3„e-*-*), 

where u = H\^ — 1/2| m. For large m, the exponentials oscillate fast compared with 
the expansion rate and we have a good definition of positive and negative frequency modes 
corresponding to our 'out' vacuum. The Bogoliubov coefficients /?„ can be readily found by 
using the expansion of the hypergeometric functions ([T2|) at late times and the normalization 
condition |anP — |/5nP = 1- The expectation value of the number of 'out' particles in a given 
mode is then 

The distribution is independent of n, as expected for nonrelativistic particles. At any given 
time, the ocupation numbers of the relativistic modes are exponentially suppressed with 
respect to the non-relativistic ones, so the distribution is cut-off at n ~ am. The resulting 
particle production rate per unit physical volume is (dN/dtdV) ~ m^H exp (— 27rm/if) ||2CI|| . 
In Ref. |19[ it was shown that the same rate can be obtained by considering the one-loop 
prefactor to the instanton contribution. There, the cut-off arises because particles with 
n >> a(t)m have not been created yet at time t. 

Hence, the tunneling boundary condition does not lead to disastrous pair production, 
but to a result which is in good agreement with the instanton calculation. 
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In following this tunneling vs. Hartle-Hawking debate, the reader should be aware that 
both wave functions are far from being rigorous mathematical objects with clearly specified 
calculational procedures. Except in the simplest models, the actual calculations of ipT and 
ipHH involve additional assumptions which appear reasonable, but are not really well justified. 
(The results presented in this paper are no exception). For a recent discussion of the problems 
associated with defining and interpreting the cosmological wave function see, e.g., Ref. |TT[ . 
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